An Integrated Chemical and Stable-Isotope Model of the Origin of Midocean Ridge Hot Spring Systems

Chemical and isotopic changes accompanying seawater-basalt interaction in axial midocean ridge hydrothermal systems are modeled with the aid of chemical equilibria and mass transfer computer programs, incorporating provision for addition and subtraction of a wide-range of reactant and product minerals, as well as cation and oxygen and hydrogen isotopic exchange equilibria. The models involve stepwise introduction of fresh basalt into progressively modified seawater at discrete temperature intervals from 100 ø to 350øC, with an overall water-rock ratio of about 0.5 being constrained by an assumed •80,2 o at 350øC of +2.0 per mil (H. Craig, personal communication, 1984). This is a realistic model because: (1) the grade of hydrothermal metamorphism increases sharply downward in the oceanic crust; (2) the water-rock ratio is high (> 50) at low temperatures and low (<0.5) at high temperatures; and (3) it allows for back-reaction of earlier-formed minerals during the course of reaction progress. The results closely match the major-element chemistry (Von Damm et al., 1985) and isotopic compositions (Craig et al., 1980) of the hydrothermal solutions presently emanating from vents at 21øN on the East Pacific Rise. The calculated solution chemistry, for example, correctly predicts complete loss of Mg and SO,• and substantial increases in Si and Fe; however, discrepancies exist in the predicted pH (5.5 versus 3.5 measured) and state of saturation of the solution with respect to greenschist-facies minerals. The calculated 6D.2 o is +2.6 per mil, in excellent agreement with analytical determinations. The calculated chemical, mineralogic, and isotopic changes in the rocks are also in good accord with observations on altered basalts dredged from midocean ridges (Humphris and Thompson, 1978; Stakes and O'Neil, 1982), as well as with data from ophiolites (Gregory and Taylor, 1981). Predicted alteration products include anfiydrite and clay minerals at low temperatures and a typical albite-epidote-chlorite-tremolite (greenschist) assemblage at 350øC. The models demand that the major portion of the water-rock interaction occur at temperatures of 300ø-350øC. Interaction at temperatures below approximately 250øC results in negative g•80.2 o shifts, contrary to the observed positive g•80 values of the fluids exiting at midocean ridge vents. Hydrogen isotope fractionation curves by Suzuoki and Epstein (1976), Lambert and Epstein (1980), and Liu and Epstein (1984), among others, are compatible with the model, and require 6D.2 o to increase at all temperatures as a result of seawater-basalt interaction.